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I. INTRODUCTION : 


L’optimisation est une branche des mathématiques. Dans la pratique, on part d’un problème concret, on le 
modélise et on le resoud mathématiquement. Des problèmes d'optimisation à objectifs multiples se 
posent dans de nombreux domaines, tels que l'ingénierie, l'économie et la logistique, lorsque des décisions 


optimales doivent être prises en présence de compromis entre deux objectifs contradictoires ou plus. 


A. PROBLEMATIQUE : 
Le secteur de la construction qui est sujet de ce rapport est pointé du doigt. Deux tiers des constructions 


de la planète sont aujourd’hui réalisés en béton armé lui-même constitué majoritairement de sable. Le 
sable dont la consommation augmente au fur et à mesure que la population et la taille des villes 
s’agrandissent, commence à se faire rare. Il devient alors crucial de trouver des alternatives au béton pour 
la construction. Face aux pénuries de matières premières, de surfaces habitables et au vu des objectifs de 
préservation de l’environnement, les techniques constructives doivent évoluer. Ces dernières nécessitent 
des usages parcimonieux de matières premières non renouvelables. Les pays industrialisés cherchent alors 
à réduire la consommation d'énergie et à utiliser des matières premières de substitution, notamment 
renouvelables. Alors que la fonction primaire d’une construction est de protéger l'usager des conditions 
climatiques extérieures et d’un environnement agressif, la conception d’un bâtiment est guidée par de 
nombreuses disciplines incluant l’énergétique, la mécanique des structures, l’environnement, l'éclairage et 
acoustique. Une construction, et donc son système constructif, doit assurer à la fois un rôle de protection, 
de sécurité et de confort de l'usager tout en limitant son impact environnemental. || s’agit ainsi d’une 
conception multidisciplinaire faisant intervenir de multiples acteurs à divers moments d’un projet de 
construction. La conception d’une construction devient un problème complexe, puisqu'il nécessite de 
traiter simultanément plusieurs objectifs parfois contradictoires. Actuellement, dans la plupart des projets, 
chacune des disciplines est traitée indépendamment des autres. Les solutions technologiques sont alors 
améliorées par recherches itératives et peu d’alternatives sont finalement évaluées. Il s’agit de tenir 
compte de compromis entre des objectifs relatifs au confort des occupants et à l'impact environnemental 
de la construction. Ces méthodes doivent si possible inclure les exigences réglementaires dès le début d’un 


projet. 
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D. SOLUTION DE LA PROBLEMATIQUE : 


Le projet de construction est un problème de conception multicritère nécessitant la prise en compte de 
performance d'usage, d'économie, de faisabilité industrielle, d'architecture et d'environnement. La 
recherche des meilleures solutions possibles constitue un problème d'optimisation multi objectif où les 
objectifs sont souvent contradictoires. L’ensemble des variables de décisions considérées définit un espace 
de solutions que le concepteur doit étudier pour chercher et trouver la meilleure solution correspondant 
aux exigences et à la volonté de la maîtrise d'ouvrage (décideur). Ce projet a pour but de définir, décrire 
et expliquer une des méthodes ď’optimisation adoptées en conception multicritère dans le domaine de 
construction et plus particulièrement en optimisation multi objectif qui est la solution à cette 


problématique. 


L Les méthodes Meta heuristiques : 


Des techniques d'optimisation robustes peuvent être utilisées lorsque dans ce domaines et se présentent 
ainsi comme solution à cette problématique. Le but est de trouver une solution réalisable pour toutes les 
données et optimale dans un sens. L'objectif est de trouver une stratégie qui soit réalisable pour toutes les 
instances de données possibles (ou presque) et d'optimiser les performances attendues du modèle. 


Les méta-heuristiques, se présentent candidat comme meilleure solution, et sont une famille d'algorithmes 
stochastiques, apparues à partir de l’année 1980, destinées à la résolution des problèmes d'optimisation 
difficiles. Elles permettent d'obtenir une valeur approchée de la solution optimale en un temps 
raisonnable. Leur particularité réside dans leur flexibilité, c’est à dire elles sont adaptables à la résolution 
d’un ensemble de problèmes dans différents domaines sans avoir à modifier le principe de base de la 
méthode (d’où le qualificatif ‘méta’). Les méta-heuristiques sont souvent inspirées de processus naturels 
qui relèvent de la physique, de la biologie de l’évolution (les algorithmes génétiques) ou encore de 
l'éthologie (les algorithmes de colonies de fourmis ou l'optimisation par essaim particulaire). L'ensemble 
des méta-heuristiques proposées dans la littérature sont partagées en deux classes : des méta heuristiques 
à base de solution unique et des méta-heuristiques à base de population de solutions. 


Les méthodes à base de population de solutions auxquelles on s’intéressera dans notre rapport : 
manipulent une population de points, et à chaque itération un nouveau jeu de points est choisi. Elles 
débutent la recherche avec une multitude de solutions. Elles s'appliquent sur un ensemble de solutions 
afin d’en extraire la plus adéquate (l’optimum global) qui représentera la solution du problème traité. 
L'utilisation d’un ensemble de solutions au lieu d’une seule solution permet une bonne diversité de la 
recherche et augmente la possibilité d'émergence de solutions de meilleure qualité. Une grande variété de 
méthodes basées sur une population de solutions a été proposée dans la littérature, commençant par les 
algorithmes évolutionnaires, passant par les algorithmes génétiques et arrivant aux algorithmes à base 
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d'intelligence par essaims (l'algorithme d'optimisation par essaim de particules, l'algorithme de colonies de 
fourmis, l'algorithme de colonies d’abeilles,.… et j'en passe. 


e CHOIX DE LA METHODE AVEC JUSTIFICATION : 


Ce projet vise à optimiser un problème de construction, y compris tous les critères de conception. C’est 
pour cette raison, j'ai mis au point une optimisation multi objectif basée sur l’OEP afin de déterminer un 
ensemble optimal de solutions. Car et comme a montré la majorité des expériences numériques réalisées 
pour un exemple concret de construction d'études de cas, la technique PSO multi-objectifs était capable de 
produire des résultats chacun d’entre eux offre une solution alternative satisfaisante pour l'optimisation 
multi-objectifs des projets de construction. De plus L’optimisation par essaim de particules se caractérise 
par sa convergence rapide, ce qui peut le conduire à stagner dans des optima locaux 


Il. PRESENTATION DE L’OFP : 
À. UN PEU D'HISTORIQUE : 


L’optimisation par essaim de particules (OEP), est un algorithme évolutionnaire largement utilisé dans les 
problèmes d'optimisation. C’est une méthode née en 1995 aux États-Unis sous le nom de Particle Swarm 
Optimization (PSO). Initialement, ses deux concepteurs, Russel Eberhart et James Kennedy, cherchaient à 
modéliser des interactions sociales entre des « agents » devant atteindre un objectif donné dans un espace 
de recherche commun, chaque agent ayant une certaine capacité de mémorisation et de traitement de 
l'information. La règle de base était qu'il ne devait y avoir aucun chef d'orchestre, ni même aucune 
connaissance par les agents de l'ensemble des informations, seulement des connaissances locales. Un 
modèle simple fut alors élaboré. Dès les premières simulations, le comportement collectif de ces agents 
évoquait celui d'un essaim d'êtres vivants convergeant parfois en plusieurs sous-essaims vers des sites 
intéressants. Ce comportement se retrouve dans bien d'autres modèles, explicitement inspirés des 
systèmes naturels. Ici, la métaphore la plus pertinente est probablement celle de l'essaim d'abeilles, 
particulièrement du fait qu'une abeille ayant trouvé un site prometteur sait en informer certaines de ses 
consœurs et que celles-ci vont tenir compte de cette information pour leur prochain déplacement. 
Finalement, le modèle s'est révélé être trop simple pour vraiment simuler un comportement social, mais 
par contre très efficace en tant qu'outil d'optimisation. Comme nous allons le voir, le fonctionnement de 
l'OEP fait qu'elle peut être rangée dans les méthodes itératives (on approche peu à peu de la solution) et 
stochastiques (on fait appel au hasard). Sous ce terme un peu technique, on retrouve un comportement 
qui est aussi vieux que la vie elle-même : améliorer sa situation en se déplaçant partiellement au hasard et 
partiellement selon des règles prédéfinies. 
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Mécanismes de communication entre chaque particule d'un essaim faisant intervenir les relations entre la mémoire 
des particules et leur position, d'après {Clerc 2014) 


B. CONCEPT GENERAL : 


+ Un algorithme PSO imite le comportement 
social d'animaux et d'insectes : 
— Les membres de la population interagissent 
entre eux tout en apprenant de leur propre 
expérience 
- Ils bougent graduellement dans des 





meilleures régions de l’espace des solutions 


+ PSO fait évoluer un essaim de particules : 
- Chacun réside à un endroit dans l’espace de recherche 


— la valeur de coût (ou fitness) lié avec chaque particule indique la qualité de sa position dans l’espace 
+ Les particules planent sur l’espace de recherche avec une certaine vélocité (direction et vitesse) 
qui est influencée par : 


1. la vitesse et direction actuelle 
2. la meilleure position qu’elle a trouvé à date 


3. la meilleure solution trouvée à date par ses voisins 
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+ L’essaim converge éventuellement à des positions ‘optimales’ 


C. EXPLICATION DE L'ALGORITHME PSO STANDARD : 





Algorithme PSO Standard (gbest) 


(1) début 
(2) pour chaque particule i 
(3) initialiser position et vélocité 


(4)  fin-pour 
(5) Tant que (non nombre maximum d'itérations atteint) 


(6) pour chaque particule i 

(7) déterminer valeur de la fonction objective y(t) 

(8) Si w(t) meilleure que w(bit)) 

(9) alors w(bi(t)) = y, et bift) = position courante 
(10) Si y(t) meilleure que w(b,(t)) 

(11) Alors y(bg(t)) = y(t) et bg(t) = position courante 
(12) fin-pour 

(13) pour chaque particule i 

(14) calculer la vélocité selon équation (1) 

(15) mettre à jour la position de la particule selon l'équation 
(3) 

(16) fin-pour 

(17) mettre à jour le facteur d'inertie selon l'équation (2) 

(18)  fin-tant que 

(19) fin 


L’algorithme du PSO est simple et se décompose en six étapes : 


1. Initialiser aléatoirement la position et la vitesse de chaque particule de l'essaim (u particules) dans 
l'espace des solutions de dimension n. 


2. Évaluer la valeur de la fonction-objectif à optimiser dans les n dimensions pour chacune des u particules 
(sa position). 


3. Comparer la valeur de chaque particule avec la meilleure valeur qu’elle n'ait jamais obtenue (pbest). Si la 
valeur actuelle est meilleure que pbest alors elle devient pbest. 


4. Identifier la particule du voisinage avec la meilleure valeur ; si cette dernière est meilleure que la 
meilleure valeur déjà obtenue par une des particules de l'essaim (gbest), elle devient la meilleure position 
jamais rencontrée par le voisinage gbest. 


5, Modifier la vitesse et la position de chaque particule selon les équations de mise à jour de la vitesse et 
de la position [V.2] et [V.3|. 


6. Retourner à l'étape 2 jusqu'à ce que le critère d’arrêt soit atteint. 


La dimension de l'espace des solutions étant fonction du nombre n de variables de l'ensemble des 
fonctions-objectifs, la position et la | = vitesse de déplacement d’une 
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particule sont des vecteurs à n composants. Les vecteurs vitesse et position sont 


actualisés 


pour chaque particule i à chaque itération t selon les équations [V.2] et [V.3] 


Ca I 
Li 


= i E M à T S T d I 2 E 
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[v.2] 


w représente le poids de linertie d'une particule ; 
r et r sont des nombres aléatoires entre 0 et 1 ; 
C1 et & sont des constantes positives d'accélération ; 


P; est la meilleure position jamais visitée par la particule i, son pbest ; 


g, est la meilleure position jamais visitée par une des voisines de la 


particule i, son gbest. 


[v.3] 


i = 1,2,..u et j = 1,2,.,n avec H le nombre de particules (taille de 
essaim) et n la dimension de l'espace des solutions ; 


Ü+ 
Xij 





+i : n na ai 
V;; est la vitesse de la particule j à l'instant t+1. 


La Figure suivante permet d'illustrer ces déplacements qui sont des sommes vectorielles pondérées. 
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III. APPLICATION DE LA PSO : 


Voici un exemple très simplifié s'inspirant d'un problème réel. On considère deux stations A et B, 
émettrices d'un champ électromagnétique sinusoïdal décroissant avec la distance. On cherche le point 
entre ces deux stations où le champ total reçu à un instant donné est minimum. La figure 4 donne une idée 
des champs en tout point entre les deux stations, distantes de 2,7 km. En posant x comme étant la distance 
à la station A, on a les équations suivantes : 


charrp _A= e "{(1+sinfiQr)} (1) 
champ_B = 6 P*A1(1 + sin(10/2,7 — z|} (2) 





On cherche le point x où la somme champ A+champ_B est minimale. On remarquera que la longueur 
d'onde commune n'étant pas en rapport simple avec la distance entre les stations, la détermination 
analytique n'est pas évidente. 


EE a 
0 M 


cn 
Ci 


_ E 
T 


Station À 
Station B 
A+B 


=É 
HI 


= 
3 ! 
C 


So 
w 


m | 
m on 


O O 
EF 


PJ 


15 


Distance à À 





De plus, les méthodes classiques (gradients, par exemple) échouent du fait de l'existence de multiples 
minimums locaux. Par contre, il est très facile d'ajouter la fonction à minimiser au programme d'OEP donné 
en annexe. Le problème est de dimension 1, et l'espace de recherche est [0; 2,7]. En moins de 200 
évaluations, on trouve x = 0,4555478 et, donc, un champ total de 0,06735316. C'est déjà une très bonne 
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approximation du minimum global. Une solution un peu meilleure peut être obtenue avec plus 
d'itérations, mais le gain est très minime. Par exemple, après 20000 évaluations, on a x = 0,4553560 et un 
champ total de 0,06735263. A titre de comparaison, avec le 


solveur d'un tableur courant du marché, on ne trouve pas du tout la solution si l'on indique le même 
espace de recherche. C'est seulement si l'on donne un point de départ déjà proche (par exemple 0,4) que 
le minimum est trouvé. Ceci suggère que, intégrée à d'autres logiciels, l'OEP pourrait être un solveur 
particulièrement efficace 


IV. ELARGISSEMENT DU SUJET : 


_ PSO a valeurs mixtes _ 


A. Problème 
Travailler avec un algorithme PSO pur implique un espace de recherche continu. Or pour des problèmes du 
monde réel, l'espace de description du problème peut contenir des descripteurs à valeurs continues, 
discrètes et/ou symboliques. 


Les données mixtes sont fréquentes dans les tâches de Résolution de problème. Par exemple, lors du tri 
des patients dans un hôpital, afin de décider s'il y a lieu de les admettre ou de les traiter en ambulatoire. Le 
résultat repose généralement sur les descripteurs qui peuvent être continus (par exemple la température), 
nominaux (par exemple la présence ou l'absence de certains symptômes) et des mesures ordinales (score 
par exemple). Un autre exemple est l'analyse d'une demande de crédit dans un contexte bancaire. 
L'évaluation de la demande se base sur des informations telles que l'âge (discret), le gain annuel (continu) 
ou le statut marital (catégoriel). 


Lorsqu'on traite des données mixtes, une manière typique de procéder passe par à une étape préliminaire 
de codage qui apparient des valeurs non numériques à des nombres entiers, mais elle se heurte à deux 
problèmes connexes : 


1) Comment établir une relation d'ordre sur les données transformées ? 


2) Comment remédier au biais causé par une représentation des connaissances plate (non pondérée) qui a 
perdu l'importance relative des différentes valeurs non numériques pendant le processus de codage. 


B. Solution : 
Notre objectif est une approche permettant de traiter les données mixtes sans les séparer. Nous 
proposons de travailler avec deux espaces (figure ci-dessous) : un espace de recherche où les particules 
évoluent de manière continue selon les lois de base de la PSO (Kennedy et Eberhart, 1995), et un espace de 
description qui reflète la réalité et où les entités peuvent avoir des valeurs mixtes. Il s'ensuit que si une 
particule évolue le long des axes continus dans l'espace des particules, ses positions sont interprétées en 
termes de descripteurs (ou attributs) de natures diverses dans l'espace de description, et une cartographie 
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sémantique entre ces deux espaces est nécessaire. Elle est assurée par des mécanismes d'interprétation 
qui font correspondre les valeurs continues dans le premier espace à des valeurs appropriées dans le 
second. Selon l'attribut qui est pris en compte dans l'espace de description, l'un des mécanismes suivants 
est utilisé pour l'interprétation des positions de l'espace des particules : 


Mécanismes d'interprétation entre l'espace de recherche et 


l'espace de description. 





i. Organigramme 


Initialiser chaque particule 


Mise à jour velocité et 
position par eq. | et 3c 


+ LS i 





Non 


Organigramme de la PSO à valeurs mixtes. 
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ii. Algorithme de la PSO a valeurs mixtes : 


En termes de pseudo-code : 


Algorithme PSO Mixte-Améliorée 


(1) début 
(2) pour chaque particule i 
(3) initialiser position et vélocité 


(4)  fin-pour 
(5) Tant que {non nombre maximum d'itérations atteint) 


(6) pour chaque particule i 

(7) interpréter la position (*) 

(8) calculer la fonction objective w(t) avec les valeurs interprétées 
(*) 

(9) Si w(t) meilleure que w(bit)} 

(10) alors w(bi(t)) = w, et bit) = position courante 

(11) Si y(t) meilleure que w(b,{1)) 

(12) Alors w(b,{t)) = yit) et b,(t) = position courante 

(13) Jin-pour 

(14) pour chaque particule i 

(15) calculer la vélocité selon équation (1) 

(16) mettre à jour la position de la particule selon l'équation (3c) 
(7) 

(17) fin-pour 

(18) mettre à jour le facteur d'inertie selon l'équation (2) 

(19)  fin-tant que 

(20) fìn 


__{ Changements apportés à l'algorithme original 


V. APPLICATION AU MAROC : 
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La forte croissance de la population mondiale favorise le développement des zones urbaines dont la 


population devrait tripler d'ici 2030. Pour éviter les conséquences néfastes d’un étalement urbain telles 


que la déforestation, la suppression de milieux naturels et l'augmentation des consommations de 


carburant liées au transport a densification urbaine est à privilégier. L’une des conséquences de la 


densification urbaine est la construction de bâtiments multiétages. Pour cela une optimisation des 


ressources est nécessaire. Le Maroc a donc eu recourt a OEP qui est reconnue au Maroc aussi bien que 


dans les autres coins du monde pendant plusieurs années comme une méta heuristique stochastique 


performante dans la résolution des problèmes d'optimisation difficiles, la méthode d’optimisation par 


om — 
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essaim particulaire « PSO » présente cependant des points faibles : le temps de calcul 
considérable et la convergence prématurée. Plusieurs études ont été menées au Maroc pour trouver le jeu 
de paramètres qui conduit à de bonnes performances de l'algorithme. Cependant cette méthode a 
apportée de nombreuses bénéfices au Maroc dans plusieurs domaines non seulement celui de la 
construction, L'algorithme proposé au Maroc a été testé afin d'améliorer la performance et la fiabilité des 
structures. 


VI. CONCLUSION : 


Les méta-heuristiques permettent l'absence d'hypothèses particulière sur la régularité de la fonction 
objective Les résultats obtenue par PSO sont très satisfaisant et confirment bien la validité de 
l'algorithme. 


le choix de paramètres reste l'un des problèmes de l'optimisation par particules d'essaim. 
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